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Pure and Random Qnantum Ising Chain : 

Shannon and Renyi entropies of the ground state via real space renormalization 

Cecile Monthus 

Institut de Physique Theorique, CNRS and CEA Saclay, 91191 Gif-sur-Yvette cedex, France 

The Shannon and the Renyi entropies of the ground state wavefunction in the pure and in the 
random quantum Ising chain are studied via the self-dual Fernandez-Pacheco real-space renormal¬ 
ization procedure. In particular, we analyze the critical behavior of the leading extensive term at 
the quantum phase transition : the derivative with respect to the control parameter is found to be 
logarithmically divergent in the pure case, and to display a cusp singularity in the random case. 

This cusp singularity for the random case is also derived via the Strong Disorder Renormalization 
approach. 


I. INTRODUCTION 

After its introduction in fluid dynamics in order to characterize the statistical properties of turbulence (see the book 
[H and references therein), the notion of multifractality has turned out to be relevant in many areas of physi cs ( see 
for instance in particular at critical points of disordered models like Anderson localization transitions I,® or 

in random classical spin models flll - [l^ . More recently, the wavefunctions of manybody quantum systems have been 
found to be generically multifractal, with many studies concerning the Shannon-Renyi entropies of the ground states 
of pure quantum spin models |19l - [3ll| . The understanding of Shannon-Renyi entropies of excited states in disordered 
models is also essential to characterize the Many-Body Localization transition [sH, . 

In the present paper, we analyze via real-space renormalization the multifractal properties of the ground state of 
the quantum Ising chain 


H = -Y^,hatat^,-Y^Kaf ( 1 ) 

i i 

where the ferromagnetic couplings Ji > 0 and the transverse fields hi > 0 are either uniform or random. For the 
pure chain, we compare with the previous results of Refs [l^ [13, [11-[1E [H, [13, [Mj in order to test the validity of the 
RG approach. For the random case, we are not aware of previous works concerning the multifractality of the ground 
state. 

The paper is organized as follows. In section m we recall the multifractal formalism for quantum wavefunctions. 
In section m we derive the real space renormalization rule for the Shannon-Renyi entropies. The application to the 
pure and to the random quantum Ising chain are described in sections IIVI and [3 respectively. Our conclusions are 
summarized in section ED 

II. REMINDER ON MULTIFRACTALITY OF MANYBODY WAVEFUNCTIONS 

For an Hilbert space of size M (growing exponentially with the volume), the expansion of a given wave-function 
|^/> > onto a basis of M vectors |m > 


M 

\i’ >= X '^rn\m > (2) 

m—\ 

involves M coefficients 'ipm subjected to the global normalization 

M 

1 =< '0I'0 >= X (3) 

m—1 

The multifractal formalism allows to characterize the statistical properties of these M weights IV'mP as we now recall. 
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A. Inverse Participation Ratios and Shannon-Renyi entropies 

The Inverse Participation Ratios Yq{M) define the exponents T{q) 

M 

Yq{M) = ^ cx 

^^ M->- + oo 

m—1 

The related Renyi entropies involve the generalized fractal dimensions Dq 

lny,(M) 


with the simple relation 


Sq(M) = ^ oc D(q) 

I - q M^+oo 

T{q) 


InM 


D{q) = 


9-1 


For g = 1, Eq. 0 corresponds to the standard Shannon entropy 


M 


^i(M) = -V l^^plnl^^P cx D{l)\nM 

^ M—>-+oo 


m=l 


(4) 

(5) 

( 6 ) 

(7) 


B. Multifractal spectrum /(a) 


Among the M configurations, the number J\fM{ce) of configurations m having a weight of order |'0mP oc M “ 
defines the multifractal spectrum /(a) 

Nm{oi) oc (8) 

M->-oo 

The saddle-point calculus in a of the Inverse Participation Ratios of Eq. |4] yields 

Yq{M) J da (9) 

where the saddle point value aq is determined by the condition 

f'{aq) = q ( 10 ) 

i.e. r(g) and /(a) are related via the Legendre transform 

T{q) + f{a) = qa 
T'{q) = a 

/'(«) = 9 (11) 

This yields the following parametric representation of /(a) 

-dq\uYq 


= ^"'( 9 ) = 


InM 


Em=i IV'mP'JlnM 


( 12 ) 


/K) 




In 




E“=ilV’».P^lnM 


(13) 


Let us now mention some important special values of the index q. 

For q = 0, the IPR of Eq. |4]simply measures the size of the Hilbert space Yo{M) = M corresponding to t(0) = — 1 
and D{0) = 1 and to /(ao) = 1 : this means that there exists an extensive number 0{M) of configurations having a 
weight M“““, so that 


MlnM 


( 14 ) 
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is called the typical exponent. 

For the Shannon value g = 1, the IPR of Eq. 0]is normalized to unity, so that r(l) = 0 and 

ai = = /(«,) = D{1) (15) 

is called the ’information dimension’ that characterizes the Shannon entropy of Eq. [T] 

In the limit q —>■ +oo, the IPR of Eq. 2] is dominated by the configuration having the highest weight Yq^+ao — 
IV’maicP'^ c>c i.e. a+oo represents the smallest exponent 


Q^+oo — 


In ^ij^rnax \ 

InM 


II(+oo) 


and the corresponding singularity spectrum vanishes /(a+oo) = 0 


(16) 


C. Application to the ground state of the quantum Ising chain 

Eor the quantum Ising chain of Eq. [T] containing N spins, the expansion of the wave function in the basis 

i^>= E E •• E (17) 

Si—il S2—il Sat— il 

involve M = 2^ coefficients i/'(5'i, S' 2 ,..., Sn)- 

1. Ferromagnetic limit hi = 0 

In the ferromagnetic limit where all transverse fields vanish = 0, the ground state is simply the linear combination 
of the two ferromagnetic states (both for periodic or free boundary conditions) 

1^ ^/erro_ l + ;+ i- -! + > 

V 2 

SO that among the M = 2^ coefficient 'i/m, only two are non-vanishing and equal to As a consequence the IPR 

of Eq. 0] reduces to 

y/"’'™(M = 2^) = 2^-9 (19) 

and the corresponding Renyi entropies all take the simple finite value 

Sf^To^M) = = in2 (20) 

^ 1 - g 

i.e. all generalized fractal dimensions of Eq. [S] vanish 

= 0 (21) 

2. Paramagnetic limit Ji = 0 

In the paramagnetic limit where all ferromagnetic couplings vanish = 0, the ground state is simply given by the 
tensor product of the ground state of cr’” for each spin (both for periodic or free boundary conditions) 

1 ^ ^para^ = -1 ^^ 2 ) 

i.e. all M = 2^ coefficients are equal to 2~^ so that the IPR of Eq. |4]read 

ypara^M = 2^) = 2^*'^“'^^ 


( 23 ) 
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and the Renyi entropies all have the same extensive value 

SP<^ra(^M) = = inM = iVln2 

i.e. all generalized fractal dimensions of Eq. [5] are equal to unity 


Tjpara _ i 
9 


(24) 


(25) 


3. Generic case 


Apart from the two extreme cases just discussed, one expects that the ground state wavefunction is characterized by 
a non-trivial multifractal spectrum. For the pure case where Ji = J and hi = h, exact expressions for the generalized 
dimension Dq as a function of the control parameter K = ^ exist in the limit q = ±oo |19l-l^ 1^ 1^ 


Dr=±o.{K) = i - 


1 


1 


2 27r In 2 


du In 


1± 


K — cos u 


y/l + K^- 2K 


cosu 


and thus also for q = 1/2 as a consequence of the duality relation [I9l - l2ll 




(26) 


(27) 


These functions are continuous in the control parameter K, but the derivatives are singular at the critical point 

K,= l 


dKDlZl^{K) 


1 


27rln2 


du 




K — cos u 


cosu 


= T- 


1 


r/C 


4K 


7r(ln2)(l + AT) V(1 + ^)V 2A:in2 


1 + — 2 AT cos 1 

-sgniK - 1 ) 


(28) 


where JC{m) = ^ is the elliptic integral of the first kind displaying the logarithmic singularity /C(m) ~ 

^ Y 1 —m sin^ 9 

— iln(l — m) for to —>■ 1“. As a consequence, the first term corresponds to a logarithmic divergence, whereas the 
second term corresponds to a discontinuity 




|if-i|<i^27r(ln2) 


In |A: - 1| + ^^sgn(A: - 1) 


(29) 


Other values of q, especially the Shannon value q = 1, have been studied numerically In the following 

sections, we describe how the Fernandez-Pacheco real-space renormalization allows to derive a simple approximation 
for the multifractal properties, both for the pure and for the random chain. 


III. ANALYSIS VIA REAL-SPACE RENORMALIZATION 

The pure quantum Ising chain is the basic model to study quantum phase transitions at zero-temperature [s^ : 
it is in the same universality class as the exactly solved 2D classical Ising model. The random quantum Ising chain 
has been solved by Daniel Fisher [ssj] via the asymptotically exact strong disorder renormalization procedure : the 
transition is governed by an Infinite Disorder Fixed Point and presents unconventional scaling laws with respect to 
the pure case. 

From the point of view of Block-Renormalization for quantum models, there exists a simple self-dual procedure 
introduced by Fernandez-Pacheco [H,|33, which is able to reproduce the exact critical point {J/h)c = I and the exact 
correlation length exponent = I of the pure chain [^, even if it is not able to reproduce the exact magnetic 

exponent [s^. Its application to the disordered chain also reproduces correctly the location of the critical point, the 
activated exponent ip = 1/2 and the two correlation length exponents even if it is not able to reproduce the 

properties dominated by rare events [s^ . In this section, we describe how it leads to simple renormalization rules for 
the Shannon and Renyi entropies of the ground state wavefunction. 
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A. Reminder on the Fernandez-Pacheco self-dual real space renormalization 

The goal is to replace each block of two spins \ <J 2 i) by a single renormalized spin ■ The idea of Fernandez- 

Pacheco [ 3 ^ , written here for the case of arbitrary transverse fields and arbitrary ferromagnetic couplings [s^ , is 

to decompose the Hamiltonian of Eq. [T] as the sum H = following choice 

Hintra — ^ ^ [ ^2i—lO'2i—1 

i 

Hextra — ^ ^ [ ^2iO’2i '^‘^i^2i^2i+l\ (^0) 

i 

Since Hintra is diagonal in the tr* basis for the even spins a^i, we have to diagonalize Hintra independently in 
each subspace {S 2 , Si, 821 , ■•)■ In addition, the odd spins a 2 i-i are independent of each other and submitted to the 
one-spin Hamiltonian 


H.[tl = - J 2 ^-lS 2 ^<Tl,_i (31) 

The two eigenvalues are independent of the value of S 2 i = ±1 

^2i-l — ^\j^2i-l + (^ 2 ) 

with the following corresponding eigenvectors 


|A2-_i(^2.)> =\ 


1 + 


J2i-lS2i 

/P ^72 
V "'22 ——1 


|S'2i_l — -|-1 > -I- 


\ 


1 - 


J2i-lS2i 




|<S'2i-l — —1 > 


|A+_i(520> =-^ 


1 - 


J2i-lS2i 

'\/^'2i-l~^’^2i-l 


|>S'2i-l — -|-1 > 


1 + 


J2i-lS2i 




|'S'2i-l — —1 > 


Projecting onto the ground state for the odd spins yields the renormalized Hamiltonian for the even spins 

^ ^~^2i^R(2i) ~ J^-2^R(2i-2)^R{2i) 
i 

with the renormalized transverse fields 


h2i = h2i 


h2i-l 


\j^ 2 i-l + 'Ali-l 


(33) 


(34) 


(35) 


and the renormalized couplings 


jR 

’J2i- 


— J2i-2 


J2i-1 


\J^2i-l + 'Ali-l 


(36) 


B. RG rules for the IPR and for the Shannon-Renyi entropies 


From the point of view of the IPR of Eq. 01 the real-space renormalization yields 


= E E E \i^{SuS2,...,SN)\^^ 

5 'i=ihlS 2 —il iSiv—il 

K 

= E E E \i^'^{S2,Si...,SN)f^f[u^i-i) 

5'2—il S4— il iSat— il i—1 


(37) 
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where yq{2i — 1; 2i) represents the IPR of the wavefunction X 2 i_i{S 2 i) (Eq. [SS]) of the odd spin 821-1 in the field of 
the even spin S 2 i 


yq{2i - l;2i) 


>S'2i-l=±l 


1 + 


J2i-1 




J2i-1 1 9 

2 


(38) 


Since yq{2i — l;2i) is actually independent of the value S 2 i = ±1, and only depends on the ratio K 2 i-i = ( , it is 

convenient to introduce the auxiliary function 


y,iK) = 


1 + 


K It 


1 - 


K 

VT+7?!' 


Then Eq. [37] can be rewritten as 


Si—il S2—il Sjsr—il 


where 


E E E \^^iS2,s^...,SN)\ 

S 2 —il S 4 —il Sjv=il 


2g 


(39) 


(40) 


(41) 


is the IPR of the renormalized system of the ^ even spins with the renormalized parameters given in Eos 155] and 1551 
The multiplicative RG rule of Eq. 001 translates into the following additive RG rule for the Shannon-Renyi entropies 
of Eq. [5] 


SqiN) = 


InYqiN) 

1-q 


= E 

i=l 


\nyq{K2i-i) , 


(42) 


where ^ entropy of the renormalized system of the ^ even spins with the renormalized 

parameters given in Eqs [OS] and jOH 

In the following sections, we analyze the consequences of the RG rule of Eq. 01] for the pure chain and for the 
random chain respectively. 


IV. MULTIFRACTAL SPECTRUM FOR THE PURE CASE 


A. RG rules for the couplings l36ll 


If the initial parameters are (/i, J) on the whole chain, the Fernandez-Pacheco renormalization yields after one RG 
step the renormalized parameters (Eqs 03] and 1551) 




h , = 

r 

V 


so that the ratio K = evolves according to the simple rule 


(43) 


(44) 


The disordered attractive fixed point K = 0 and the ferromagnetic attractive fixed point K —>■ -boo are separated by 
the unstable fixed point = 1 characterized by the correlation length exponent v = 1 obtained by 2v = (j)'(^K,,) = 
2Kc = 2. So the location of the critical point Kc = 1 and and the correlation length exponent 1 / = 1 are in agreement 
with the exact solution [s^. 
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B. RG rule for the Shannon-Renyi entropies 


The factors yq{K 2 i-i) in Eq. [38] are all the same, and depend only on the control parameter K = J/h so that the 
RG rule of Eq. 021 for the Shannon-Renyi entropies reads 


where the function yq{K) has been introduced in Eq. 1391 The iteration yields 

NlnyqjK) N\nyq{K^) N \nyq{K^) N \nyq^ 


Sq{K-N) = 


2 l-q 4 1-9 


1-9 


2^= 1-9 


(45) 


(46) 


where the term corresponds to the contribution of the RG step, where ^ spins coupled via are 

eliminated, and where ^ spins survive with the renormalized coupling . 


C. Generalized fractal dimensions Dq{K) 


Erom Eq. HU one obtains the generalized fractal dimensions Dq{K) of Eq. [5] as a series over the RG steps 
A: = 1, 2,.., - 1-00 


Sq{K-,N) 


Dq{K) = lim 


N^+ca A'^ln2 2^ 

k=l 


in terms of the auxiliary function dq{K) obtained from Eq. | 


dq{K) = 


\nyq{K) 


In 


i-t 






(1 — 9 ) In 2 

that becomes for the Shannon value 9 = 1 


(1 — 9 ) In 2 


dq^l{K) =-^ 


-1 1 K ■ 

2 

In 

■1 1 K ■ 

2 

-b 

K ■ 

Vi + k2 

2 

In 

K - 

2 


In 2 

As it should, one recovers at the paramagnetic fixed point K — 0 where dq{K = 0) = 1 the value of Eq. 

Dq{K = 0) = 1 

and at the ferromagnetic fixed point K = -l-oo where dq{K = -too) = 0 the value of Eq. |2I] 

Dq{K = -l-oo) = 0 


(47) 


(48) 


(49) 

(50) 

(51) 


D. Multifractal spectrum at the critical point Kc = 1 

At the critical point Kc = 1, the generalized fractal dimensions Dq{Kc = 1) of Eq. 07] simply coincide with the 
auxiliary function of Eq. 08] 


Dq{Kq = 1 ) 


= d,{l) 



(1 — 9 ) In 2 


(52) 



















































FIG. 1: Pure quantum Ising chain at the critical point Kc = 1 

(a) Generalized fractal dimension Dq{Kc = 1) as a function of the Renyi index q 

(b) Gorresponding multifractal spectrum f{a) as a function of a. 


This function is plotted on Fig. [T] (a) for the Renyi index in the interval q G [—5,+5]. Let us quote some special 
values 


Dg^o{K = l) 

Dq^l/2{K = 1 ) 


= l-| + 0(g2) 

In 2 


0.771553 


Dq=^{K=l) 

Dq=2{.K=l) 



In4-ln3 


In 2 


In 2 

0.415037 


0.600876 


Dq=+^{K=l) 


In 




In 2 


0.228447 


(53) 


in order to compare them with the explicit exact values for g = 1/2 and q = +oo (^ . and the numerical values 
for g = 1 and g = 2 [2^ 


Tiexact / 

1 ) 

2Catalan 

= , „ = 0.841267.. 

7rln2 


b 

II 1 
(/S 

II 

1 ) 

= = 0.61065... 

In 2 


b 

II 1 

II 

1 ) 

= = 0.30845... 

In 2 


r^exact / 
-^Q^+ooV-^ ~ 

1 ) 

= 1 2 ^“'“'“" = 0 . 158733 .. 

7rln2 

(54) 


Our conclusion is thus that the RG yields a very good approximated value for the Shannon value g = 1, whereas the 
approximation is clearly worse for the other Renyi indices g 1. This seems to indicate that the Fernandez-Pacheco 
RG procedure is well suited for the initial Hamiltonian corresponding to g = 1, whereas it should probably be modified 
for g 1 in order to take into account the deformation of the measure with respect to the initial quantum problem 
corresponding to g = 1. Indeed from the point of view of the e quiv alent 2D classical Ising model, the case where 
M = 2g is an integer corresponds to the ’Ising book’ (see Fig. 4 of [^) where M = 2q half-planes are glued together. 
This clearly changes the measure of the spin distribution and should be taken into account within the RG procedure 
in order to obtain better values for g I. 

On Fig. [T](b), the corresponding multifractal spectrum fKc=i{c‘) is plotted as a function of a (see section |TTB] for 
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a reminder on its meaning) from the parametric form of Eqs [12] and [T3| 

ccqiKc = 1) = — 


1 

/ 

2 

Q 

In 

2 


^ ^2 

2 

q 

In 

^ ^2 

2 

\ 

In 2 

V 

1+7^1 

2 

Q 

+ 

^“771 

2 

q 

/ 


In 


fK^=i{aq) =qaq{Kc = l) + 


■/2 


+ 


1+^ 


In 2 


In particular for <7 = 0, the typical exponent of Eq. 1141 takes the simple value 

ao{Kc = 1 ) =1 


(55) 


(56) 


E. Singularities at the critical point Kc = 1 




FIG. 2: Pure quantum Ising chain at the Shannon value q ~ 1 

(a) The generalized fractal dimension Di{K) is continuous a function of the control parameter K — 

(b) The derivative dKDi{K) as a function of the control parameter K — displays a logarithmic singularity at the critical 
point Kc = 1. 


The values Dq{K = 1) are not universal, but the derivatives of Dq{K) with respect to the control parameter K are 
expected to display universal singularities at the critical point. The first derivative of Eq. HTl 


1 +00 

i X -^ , 0^-1 




(57) 


k=l 


clearly diverges a,i Kc = 1. Since this divergence is due to the region of large RG-step k, we may evaluate the form 
of the singularity by replacing the series over the integer k by the corresponding integral over the real variable k, and 
then perform a change of variable from k to Kji = 


D'{K) 


|tf-i|«i 2 

1 f 




dkK^ d'JK^ ) 


+ 00 


dKn 


-d'JKu) 


|if-i|«i2y^ (lnif^j)(ln2) 
i.e. the second derivative displays a pole singularity 

d'qjK) <(1) 


(58) 


D''{K) 


|if-i|«i 2(ln2)(lnR:) |/f-i|«i 2(ln2)(R:-l) 


+ ... 


(59) 
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The singularity of the first derivative is thus given by the following logarithmic divergence 

<( 1 ) 


DM 


|x-i|<i 2(ln2) 


ln\K- 11 + ... 


(60) 


This behavior is in agreement with the logarithmic singularity of Eq. [2^ concerning the limit q —>■ ±oo. 
Finally, the function itself is continuous with the following logarithmic singularity 


<(1) 

DJK) ~ DM) - -P—- 
^ |K-i|«i ^ 2 ln2 


{K-l)ln\K- 1| + ... 


(61) 


For the Shannon index q = 1, the generalized fractal dimension Di{K) and its derivative dxDi^K) are plotted as 
a function of the control parameter K — on Fig. [2J The corresponding coefficient of the logarithmic singularity on 
these figures is 


arcco 


ith(-\/2) 


2(ln2) 4\/2(ln2)2 


= -0.32429 


(62) 


F. Duality functional relation 

As a final remark, let us mention that the RG result of Eq. 07] satisfies the exact duality relation quoted in Eq. \T7\ 
Indeed for the index value q = 1/2, the auxiliary function of Eq. |48| reads 


In 


dn^hiK) = 


1 + 




1 - 


y/l + K"^ 


/ 1 [ M 
In 


ln(l + 7T^) 


In 2 


In 2 


1 - 




+ 2 


1 - 


l + K^ 


In 2 


(63) 


whereas for the index value q = +oo, it reads 


In ( 


dq—-i^oo (AT) — 


In 2 


^ "^9=^ ( K 


(64) 


so that the duality relation of Eq. [571 is satisfied term by term in the series representation of Eq. 07] 


V. MULTIFRACTAL SPECTRUM FOR THE RANDOM CHAIN 
A. RG rules for the couplings |39l . l4Cll| 


In terms of the ratios 


K^ = 


Ji 

hi+i 


(65) 


the RG rules of Eq. 05] and 06] read 


K 


R 

2i-2 


jR 

_ ^2i-2 

“ hR 

"'2i 


J2i-2J2i-l 

^2i-1^2i 


K2i-2K2i-l 


that generalize the pure RG rules of Eq. 03] 

Since this corresponds to a simple addition in log-variables 


( 66 ) 


In Ar^_2 = In K2^-2 + In K2i-i 


( 67 ) 
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one obtains after k RG steps in terms of the initial variables 

InRT'^ = ^^InRTi-i = [In Ji_i — In/i^] 

i=l i=l 


The Central Limit theorem yields the asymptotic behavior 

InRT'^ ~ 2^ lln Ji_i — In/id + V^-\/[Rar [In Jd + Var\\nhiV[u 

k—y+oo ^ ■' 


( 68 ) 


(69) 


where u is a Gaussian random variable of zero mean and of variance unity. So for large RG step fc, the probability 
distribution Pk{Kii) of the renormalized coupling is log-normal 


Pk{KR)^ ~ 


(In In Ky 


k^+oo Kii'J2ttV2^ 


(70) 


with the notations 


IniL = In Ji_i — In hi 

V = Far [In Ji] -|- Far [In hi] 

The first term of Eq. [69] yields that the critical point corresponds to the condition 


In Kc = In Ji_i — In hi =0 

and that the typical correlation length exponent with respect to the length scale L = 2^ is 


(71) 


(72) 


^typ — 1 ('^^) 

Outside criticality, the competition between the first and the second term of Eq. [69] shows that the finite-size 
correlation exponent is 


ly^S = 2 (74) 

At criticality where the first term vanishes, the second random term of order corresponds to an Infinite Disorder 
Fixed Point of exponent 


^ ^ (75) 

All these conclusions of Eqs |72l IT^ [75] obtained via the application of the Fernandez-Pacheco renormalization 
to the random quantum Ising chain [39l l40| . are in agreement with the Fisher Strong Disorder renormalization exact 
results [ 3 ^ . 


B. RG rule for averaged Shannon-Renyi entropies 


Averaging Eq. [42] over the disorder yields the following RG rule for averaged Shannon-Renyi entropies 


Sy{N) 


N\nyq(K2^-l) , 
~2 + 


(76) 


where S^{^) is the averaged entropy of the renormalized system of the ^ even spins with the renormalized parameters 
given in Eq 1661 Note that whenever multifractality appears in disordered systems (see for instance the review 0 
on Anderson localization transitions), one should make the distinction between the ’averaged multifractal sprectrum’ 
based on the disorder-average of the IPR, and the ’typical multifractal sprectrum’ based on the disorder-average of 
the logarithm of the IPR. Here we only consider the disorder-averaged entropies of Eq. [TSl corresponding to the 
disorder-average of the logarithm of the IPR that define the typical multifractal spectrum. 
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FIG. 3: Random quantum Ising chain at the Renyi index q = 1 

(a) Information dimension Di{K) as a function of the control parameter K = 

(b) The derivative as a function of the control parameter K displays a cusp singularity at the critical point Kc = 1. 


C. Generalized fractal dimensions D„ 


So in terms of the auxiliary function dq{K) of Eq. 


Da 


= lim 




the iteration of Eq. [76l yields the generalized fractal dimension 

dq{Kli) 


N—y+oo In 2 


+ 00 

E 

k^l 


r*+oo 


dKiiPk-i{Kii)- 


2 k 


(77) 


where PkiKn) is the probability distribution of the renormalized variable Kr after k RG steps (Eq. [S51) . Of course 
the values Dq of Eq. [^are non-universal as in the pure case, since they depend on the contribution of the first RG 
steps. In particular they thus depend upon the initial disorder distribution Po{K). On Fig. [U we show Di{K) and its 
derivative OrDi^K) as a function of the control parameter K of Eq. [HI for the special case where the initial disorder 
distribution Pq{K) is itself given by the log-normal form of Eq. (TO] with the value V = 1, so that Eq. [TDlis valid for 
all RG steps k 


n + OO 


Dq{K) = 


+ 00 

E 

+ 00 /* + oo 


dKn 


dKn _ (lnKH-2>=-MnK)2 dq{K r) 

-2V2fc-l —1--- 


KrVQjtVW^^ 


2 k 


= E 


du _ 2 iLdq{K‘^'" ^ 
■.e 2 




2 k 


(78) 


On Fig. [21 it is clear that Di{K) and its derivative dKDi{K) are both continuous at the critical point = 1, but 
that the derivative displays a cusp (see the difference with the corresponding Figure [2] concerning the pure case). To 
understand the origin of this singularity, we may take into account the Infinite Disorder nature of the critical point 
and the Strong Disorder nature of the critical region. More precisely at large RG step k, the renormalized coupling 




(79) 


is either very small (if 2 ^ ^\nK + WoJ^u < 0) or either very big (if 2^ ^\nK + WoJ^u > 0) with the following 
leading contributions of the auxiliary function of Eq. 05] 


dq{K) ~ 1 + 0(K^) 

^ K^O 

Let us thus make the approximation that the singular part is given in the critical region by 


+ 00 




k^l 


2 k 


V du _ 1 

- 1 =^ 

o V^ 


(80) 


(81) 
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Since the singularity is due to the region of large RG-step A:, we may evaluate the first derivative by replacing the 
series over the integer k by the corresponding integral over the real variable k, and then perform a change of variable 
from k to X = ^^^(Ini^)^ 




+ 00 




1 


/• + 00 


~ -2 


V2ttVW^K 
1 




dk 


I In if I 


V2ttV2>^-^K 
/•+00 


-(JnKf 


dx 

{hv2)Vy/^K l <.lnK)2 

1 


(ln2)yv^if 


2VW-2y/^\ In if I + 2 ^^^^ + O ((lnif)2) 


(82) 


is continuous at the critical point ifc = 1, but displays a cusp as a consequence of the second term containing the 
factor I In if I oc |if — 1| near criticality. 


D. Analysis of the Shannon-Renyi entropies via Strong Disorder RG 


For the random quantum Isin g ch ain, the strong disorder renormalization procedure has been shown to be asymp¬ 
totically exact by Daniel Fisher |35| so that it can be used to derive all universal critical properties (see and the 
review [11), as well entanglement properties (see [l^ and the review (i^). It is thus interesting to consider its ap¬ 
plication to the Shannon-Renyi entropies, even if only the singularities will be meaningful, whereas the non-universal 
values dominated by the first RG steps will of course not be reproduced faithfully. 

In the Strong Disorder RG procedure [H, |4l| , the strongest term V, among the ferromagnetic couplings Ji and the 
transverse fields hi 


n = max(Ji, hi) 


(83) 


is decimated iteratively. If the strongest term is a coupling Jj, the two spins linked by this coupling are coupled 

ferromagnetically to form a new renormalized spin (the two kept states are | -I- -f > and |->). So the contribution 

of this decimation to the Shannon-Renyi entropies is actually zero 

5"“*™(iV) = S^{N - 1) (84) 

If the strongest term is a transverse field hi^ the corresponding spin is in the eigenstate of cr®, so that the contribution 
of this decimation to the Shannon-Renyi entropies is the maximal value In 2 

5'““™ (TV) = S^{N - 1) -k (In 2) (85) 

As a consequence, the generalized fractal dimension is simply 


, SgjN) Mh 

N\n2 Nh+Mj 


( 86 ) 


where J\fh is the total number of /i-decimations along the RG flow, and A/j is the total number of J-decimations along 
the RG flow. 

We refer to Refs |4lj for the renormalization rules and for the asymptotic probability distributions of the 
renormalized couplings and of the renormalized transverse fields as a function of the RG scale F = — In fl and of the 
usual control parameter 5 [H, HTj related to the previous notations {K, V) of Eq. [7T] 


S = 


InJi-i—lnhi In AT 

Var\biJi] + Var\hYhi] V 


(87) 


Here we only need the probability Fo(r)dr to decimate a coupling in the interval [F,r -|- dF] and the probability 
i?o(r)dr to decimate a transverse field in the interval [r,F + dF], that read [3^.1^ 


^o(r) = 
i?o(r) = 


2S 
1 — 

2S 

g2Sr _ I 


( 88 ) 
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In particular, the density n 5 (r) of surviving spins at scale F evolves according to 

arn5(r) = - [Po(r) + i?o(r)] n5(r) 

and the solution starting from the initial condition n 5 (ro) = 1 at the initial scale Fq reads 

t-r\ _ sinh^^Fo 
“ sinh^ ST 

The fraction of /i-decimations is then given in the thermodynamic limit by 

'To 




Mh + A/j 
= sinh^ (5Fi 


1 - 25 r 

2^ 


- 1 


sinh^ (5F 


r=+oo 


r=ro 


At criticality (5 = 0, the fraction of /i-decimation is of course 


(89) 

(90) 


(91) 


= =- 


(92) 


by symmetry. In the ordered phase (5 > 0, it is given by 

Dq{5 >0) = 1 - (93) 

whereas in the disordered phase 5 < 0 it reads 

Dq{S < 0) = (94) 

In particular the derivative 

dsDqiS) = Foe-2|^|ro (95) 

is continuous but displays a cusp at criticality 5 = 0, in agreement with the result |82] of the Fernandez-Pacheco 

renormalization. 


VI. CONCLUSION 

In this paper, we have studied via real space renormalization the Shannon and the Renyi entropies of the ground 
state wavefunction in the pure and in the random quantum Ising chain. Our main conclusion is that the leading 
extensive term presents the following singularities at the quantum phase transition : the derivative with respect to 
the control parameter is logarithmically divergent in the pure case, and displays a cusp singularity in the random case. 
This cusp singularity for the random case has been also confirmed via the Strong Disorder Renormalization approach. 
Although the logarithmic divergence for the pure quantum chain is reminiscent of the logarithmic divergence of the 
specific heat of the 2D classical Ising model, the cusp singularity found for the random quantum chain is different 
from the specific heat singularity of the corresponding 2D McCoy-Wu model (4l| . 

We should stress that the leading extensive terms of the Shannon and the Renyi entropies do not depend on the 
boundary conditions. As other extensive thermodynamic observables, their values are dominated by the small scales 
in the bulk and are thus non-universal, whereas their singularities at the critical point are dominated by the large 
scales and are thus universal. The subleading non-extensive terms of the Shannon and the Renyi entropies for finite 
chains of size N are of course also very interesting. In the pure chain, these subleading contributions are universal, 
but depend upon the boundary conditions as follows : 

(i) for the periodic chain, the subleading term is of order 0 ( 1 ) and is non-trivial only at the critical point = 1 

and for the Shannon value q = l [l9l - l^ . [241^ . , whereas the direct application of the real space RG described in 

the text would produce a non-trivial 0(1) term only at the critical point = 1, but for all value of the Renyi index 
q. This suggests, as already mentioned after Eq. that the RG procedure probably needs to be improved for (7 ^ 1, 
to take into account the deformation of the measure with respect to the initial quantum problem corresponding to 
q=l. 

(ii) for the open chain, the subleading term at criticality is logarithmic and q-dependent [ 2 ^ [H, js^, HH , as predicted 
from the presence of ’corners’ in the Conformal Field Theory of the replicated associated statistical physics model. 
It would be interesting to better understand the origin of these logarithmic terms directly in the framework of the 
quantum Ising chain, and to reproduce them by some appropriate real-space RG procedure. 
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